Relativistic many-body calculations of electric-dipole matrix elements, lifetimes and 

polarizabilities in rubidium. 
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Electric-dipole matrix elements for ns — n'p, nd — n'p, and 6d — 4/ transitions in Rb are calculated 
using a relativistic all-order method. A third-order calculation is also carried out for these matrix 
elements to evaluate the importance of the high-order many-body perturbation theory contributions. 
The all-order matrix elements are used to evaluate lifetimes of ns and np levels with n = 6, 7, 8 and 
nd levels with n = 4, 5, 6 for comparison with experiment and to provide benchmark values for 
these lifetimes. The dynamic polarizabilities are calculated for ns states of rubidium. The resulting 
lifetime and polarizability values are compared with available theory and experiment. 
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I. INTRODUCTION 

A recent proposal jlj] for quantum computation utilizes 
the ground hyperfine states of a neutral atom as qubits, 
and realizes two-qubit quantum gates by conditional ex- 
citations to Rydberg states. The atoms are confined to 
the sites of an optical lattice, and the optical potential 
seen by the atom depends on the dynamic polarizability 
of the atom. Therefore, the atom will generally move in 
a different lattice potential when it is excited to the Ry- 
dberg state during the gate operation, which may cause 
motional heating and lead to decoherence. In a recent 
paper 1^, we proposed two solutions to this problem by 
matching the ac polarizabilities of the atom in the ground 
and Rydberg states. In the first scheme, the polarizabil- 
ities are matched for the specific values of the lattice 
photon frequency between the resonances. In the sec- 
ond scheme, some accidental matches between transition 
energies are used to match the ground state polarizabil- 
ity with the polarizability of the selected Rydberg states 
(lip and Ibp in the case of Rb). The elimination of mo- 
tional decoherence is important in helping to design a 
high-fidelity two-qubit gate capable of meeting the error 
threshold for scalable quantum computation. 

Despite the existence of high-precision measurements 
of the primary transition electric-dipole matrix elements 
in alkali- metal atoms, accurate experimental data for 
other transitions are lacking with the exception of a very 
few transitions in Cs owing to the study of parity non- 
conservation. In light of the importance of the atomic 
calculations for the quantum logic gate scheme with con- 
ditional excitations to Rydberg levels, we have calculated 
the electric-dipole matrix elements for subsequent evalu- 
ation of lifetimes and polarizabilities for a number of Rb 
levels. We note that the interest in this particular logic 
gate scheme with neutral atoms results from its potential 
for fast (sub-microsecond) gate operations. 

We have performed all-order calculations of the 7s— np, 
8s— np, Ad—np, 5d—np, and 6(i—np electric-dipole matrix 
elements with n = 5, 6, 7 and 8, and of the 6d — 4/ ma- 
trix elements. These results are combined with previous 



all-order calculations of the 5s — np and 6s — np matrix 
elements Q to obtain values of the lifetimes of the 6s, 
7s, 8s, 6p, 7p, 8p, Ad, 5d, and 6d levels. The lifetime of 
the 6(^3/2 level is of special interest owing to the large 
discrepancies between existing experiments 0, The 
third-order matrix-element calculation has also been car- 
ried out to evaluate the importance of the higher-order 
contributions. The dynamic polarizabilities of the Rb ns 
states are also calculated for both low-lying and Rydberg 
levels. The evaluation of the accuracy of the polariz- 
ability calculations is conducted including a comparison 
between our data and other theory and experiment. 



II. ELECTRIC-DIPOLE MATRIX ELEMENTS 

The electric-dipole matrix elements for the 5s— bp tran- 
sitions were measured to high accuracy in |^ and ma- 
trix elements for ns — n'p transitions with n = 5,6 and 
n' = 5, 6, 7, 8 were calculated using a single-double (SD) 
all-order method in Refs. 0j[3- I'^ the present work, we 
have calculated matrix elements for 7s — np, 8s — np. 
Ad — np, 5d — np, and 6d — np transitions with n = 5, 6, 7 
and 8 as well as 4/ — 6d transitions using SD all-order 
method Q . Very briefly, the wave function is represented 
as an expansion 
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where <I>^ is the lowest-order atomic state function, which 
is taken to be the frozen-core Dirac-Hartree-Fock (DHF) 
wave function of a state v; aj and a; are creation and 
annihilation operators, respectively. The index a is 
used to represent core states and indices m and n indicate 
excited states. 

The quantities pma and Pmv are single core and va- 
lence excitation coefficients and the quantities Pmnab and 
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TABLE I: Absolute values of all-order SD electric-dipole matrix elements for the np — n's transitions in Rb in atomic units. 
The lowest-order (DHF) and third-order values all also given. All-order values for 5s — np and 6s — np transitions are from 
Ref. 0. The corresponding reduced oscillator strengths are given in rows labeled /. All-order matrix elements and experimental 
energies are used in calculation of oscillator strengths. 





5s - 5pi/2 


5s - 6pi/2 


5s - 7pi/2 


5s - 8pi/2 


5s - 5p3/2 


5s - 6p3/2 


5s - 7p3/2 


5s - 8p3/2 


Drit 


4.819 


0.383 


0.142 


0.078 


6.807 


0.606 


0.237 


0.136 


inird order 


4.181 


0.363 


0.130 


0.069 


5.899 


0.583 


0.224 


0.125 


A 11 J ,,,, 

Ail order 


4.221 


0.333 


0.115 


0.059 


5.956 


0.541 


0.202 


0.111 


J 5s^np 


U.04U4 


A A A A A 
0.004U 


A AAAC; 

O.OOOd 


A AAAO 
0.0002 


A CAAC 

0.6900 


A A1 na 
0.0106 


A AA1 '7 
O.OOi ( 


A AAA^; 

0.0006 




6s - 5pi/2 


6s - 6pi/2 


6s - 7pi/2 


6s - 8pi/2 


6s - 5p3/2 


6s - 6p3/2 


6S - 7p3/2 


6s - 8p3/2 


Urit 


4.2o6 


1 A nod 
lO.ZOD 


A ATC 

0.9(6 


A OTC 

0.3(0 


6.18( 


1/1 AK^ 

14.45 ( 


1 /I AO 


A CAT 

0.09( 


iliird order 


4.189 


9.584 


1.050 


0.420 


6.115 


13.447 


1.610 


0.668 


All order 


4.119 


9.684 


0.999 


0.393 


6.013 


13.592 


1.540 


0.628 


J6s^np 


-U.i946 


A C 1 A1 

O.olOl 


A A1 1 
0.011 / 


A AAOO 


A Am n 
-0.4Ui9 


1 Aocy 
1.026/ 


A AO'yA 
0.02/9 


A AAC O 

O.OOoo 




7s - 5pi/2 


7s - 6pi/2 


7s - 7pi/2 


7s - 8pi/2 


7S - 5p3/2 


IS - 6p3/2 


7S - 7p3/2 


7s - 8p3/2 


DHF 


0.981 


9.360 


17.612 


1.801 


1.393 


13.552 


24.708 


2.728 


Third order 


0.952 


9.304 


16.679 


1.944 


1.347 


13.517 


23.349 


2.943 


All order 


0.954 


9.189 


16.844 


1.865 


1.352 


13.353 


23.587 


2.833 




-0.0190 


-0.3330 


0.6565 


0.0186 


-0.0375 


-0.6821 


1.3170 


0.0432 




8s - 5pi/2 


8s - 6pi/2 


8s - 7pi/2 


8s - 8pi/2 


8s - 5p3/2 


8s - 6p3/2 


8s - 7p3/2 


8s - 8p3/2 


DHF 


0.514 


1.922 


16.151 


26.817 


0.727 


2.705 


23.343 


37.577 


Third order 


0.500 


1.839 


16.162 


25.587 


0.705 


2.578 


23.428 


35.770 


All order 


0.504 


1.853 


15.982 


25.831 


0.710 


2.600 


23.171 


36.123 


fss—^np 


-0.0063 


-0.0278 


-0.4701 


0.7987 


-0.0124 


-0.0539 


-0.9595 


1.5993 



Pmnva are double core and valence excitation coefficients, 
respectively. We obtain the equations for the excitation 
coefficients by substituting the above wave function into 
the many-body Schrodinger equation. The equations are 
solved iteratively and the resulting expansion coefficients 
are used to calculate matrix elements. In such a proce- 
dure, certain classes of MBPT terms are summed to all 
orders. However, the restriction of the expansion to 
the single and double excitations leads to some missing 
terms in the expression for the matrix elements start- 
ing from the fourth order. As it is important to under- 
stand how significant higher-order corrections are for the 
transitions under consideration, we also conduct separate 
third-order calculation of these matrix elements (includ- 
ing random-phase approximation (RPA) contributions it- 
erated to all orders). The third-order calculation follows 
that of Ref. ■ Such a calculation is of lower accuracy 
than the all-order one but the difference in third-order 
and all-order values provides an estimate of the impor- 
tance of the higher-order contributions. Also, the break- 
down of both all-order and third-order calculations to 
different contributions yields information regarding rela- 
tive importance of the specific terms and possible large 
cancellations between different terms. Such information 
can be used to approximate some of the omitted contri- 
butions and further estimate the uncertainty of the all- 
order calculation. We should note that despite the fact 
that the entire third-order contribution is contained in 
the all-order result, the extraction of the third-order part 
from the all-order calculation is not a straightforward 
task (see Ref. jJl fo^' the correspondence of the terms), 
so a separate calculation is made to obtain third-order 



values. The separate calculation also allows us to iterate 
RPA contribution to all orders. 

We use the system of atomic units, a.u., in which e, me, 
47reo and the reduced Planck constant h have the numer- 
ical value 1. Polarizability in a.u. has the dimensions 
of volume, and its numerical values presented here are 
thus measured in units of Og, where oq ~ 0.052918 nm is 
Bohr radius. The atomic units for a can be be converted 
to SI units via a/h [Hz/(V/m)2]=2.48832x lO'^a [a.u.], 
where the conversion coefficient is ATreoag/h and Planck 
constant h is factored out. The atomic unit of frequency 
Lu is Eh/h ~ 4.1341 X 10^^ Hz, where Eh is Hartree en- 
ergy. 

The third-order and all-order results for ns — n'p tran- 
sitions with n, n' = 5, 6, 7 and 8 are summarized in 
Table together with the lowest-order Dirac-Hartree- 
Fock (DHF) values. The corresponding reduced oscillator 
strengths calculated using formula 

fa^ = ~7l^^MD\\b)\' (2) 
[2ja + 1)A 

are listed in the rows labeled fns—>n'p- In Eq- 0, D 
is the dipole operator and A is a transition wavelength 
in nm. The all-order matrix elements and experimental 
energies from are used in the oscillator strength cal- 
culation. The sum of the f^s^np oscillator strengths with 
n = 5,6,7,8 slightly exceeds one because the contribu- 
tions to the sum rule fss^np = 1 from the transitions 
with n = 2, 3, 4 are negative. 

The results for the nd — n'p matrix elements with 
n = 4, 5, 6 and the n' = 5, 6, 7, 8 and 6d — 4/ matrix ele- 
ments are summarized in Table IhI The transitions con- 
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TABLE II: Absolute values of all-order SD electric-dipole 
matrix elements for the np — n'd and 6d — 4/ transitions in 
Rb in atomic units. The lowest-order and third-order values 
all also given. In the cases where lowest order and SD data 



are of different 


sign relative 


signs are 


shown. 




n 


5 


6 


7 


8 






4^3/2 - 


- npi/2 




DHF 


9.046 


6.725 


1.181 


0.578 


Third order 


8.092 


5.289 


1.119 


0.566 


All order 


7.847 


4.717 


1.054 


0.541 






dd'i /o - 






DHF 


4 (189 


Zr. you 


U.OO'4: 


0.262 


Third order 


3.655 


2.307 


0.502 


0.255 


All order 


3.540 


2.055 


0.470 


0.242 






4d^ It - 








1 9 941 


8 S9Q 


± .uu± 


n 787 

u. * * 


Third order 


10.964 


6.915 


1.505 


0.764 


All order 


10.634 


6.184 


1.411 


0.726 






5(^3/2 ~ 






DHF 


0.244 


18.701 


13.639 


2.660 


Third order 


1.220 


18.241 


10.600 


2.573 


All order 


1.616 


18.106 


9.768 


2.400 






5(^3/2 - 


- "■P3/2 




DHF 


0.157 


8.443 


5.983 


1.198 


Third order 


0.607 


8.232 


4.609 


1.146 


All order 


0.787 


8.160 


4.242 


1.064 






5^5/2 - 


- np3/2 




DHF 


0.493 


25.340 


17.884 


3.592 


Third order 


1.821 


24.695 


13.843 


3.435 


All order 


2.334 


24.491 


12.798 


3.201 






6^3/2 " 


- npi/2 




DHF 


0.512 


-0.254 


31.349 


22.584 


Third order 


1.076 


1.375 


31.563 


17.647 


All order 


1.180 


1.989 


31.422 


16.631 






6^3/2 - 


- 'T-P3/2 




DHF 


0.255 


-0.028 


14.158 


9.901 


Third order 


0.513 


0.732 


14.242 


7.662 


All order 


0.558 


1.012 


14.161 


7.215 






6(^3/2 - 


" np-i/2 




DHF 


0.778 


-0.047 


42.500 


29.603 


Third order 


1.532 


2.184 


42.709 


23.041 


All order 


1.658 


2.974 


42.481 


21.784 


6q 


's/2 - 4/5/2 6^5/2 - 4/5/2 6^5/2 - 4/7/2 




DHF 


6.109 


1.642 


7.343 




Third order 


9.150 


2.443 


10.924 




All order 


9.938 


2.642 


11.813 





taining the fine structure components, such as ns — 5pi/2 
and ns — 5p3/2 pairs, for example, have very similar rela- 
tive correlation contributions so we will omit the angular 
momentum j subscript in the subsequent discussion. We 
find that the relative correlation correction contribution 
varies very significantly with the transition. It is small, 
below 7%, for all 6s — np, 7s — np and 8s — np transitions 
and for some nd — n'p transitions. In those situations 
where the third order and all-order values differ by less 
than a few percent we expect the dipole matrix elements 



to be accurate to at least 2%. 

We should note that in some cases it is possible that 
such good agreement of the third-order and all-order val- 
ues is fortuitous and may result from accurate cancel- 
lations of the high-order terms. We address this issue 
below in more detail. 

In Ref. 0, the all-order values for the primary tran- 
sitions were found to agree with a recent high-accuracy 
measurement to within 0.2%-0.35%. The accuracy 
of the 5s — 6p, 5s — 7p, and 5s — 8p matrix elements 
is, however, substantially lower since they are relatively 
small matrix elements with significant relative correla- 
tion contributions. To investigate the correlation con- 
tributions in more detail, we give a breakdown of both 
third-order and all-order calculations for eight selected 
transitions in Table IIIII The third-order contributions 
are separated to random-phase-approximation (RPA), 
Brueckner-orbital (BO), structural radiation (SR) and 
normalization (Norm) corrections following Ref. For 
the all-order contributions, three terms (a, c and d) are 
listed separately and the other 17 terms are grouped to- 
gether in row labeled "Other" . The derivation and ex- 
pressions for these terms are given in Refs. 0, A 
normalization correction is given in row labeled "Norm" . 
As it was noted earlier, there is no straightforward cor- 
respondence of the all-order and third-order breakdown 
(full description is given in 0), but term a partly cor- 
responds to RPA-like corrections and term c to BO-like 
corrections. Term d is normally small with exception of 
some transitions. It is quadratic in single-valence exci- 
tation coefficients pmv and, therefore, contains only fifth 
and higher-order terms. We find strong cancellations be- 
tween BO and RPA terms for 5s — 6p, 5s — 7p, and 5s — 8p 
matrix elements. The difficulty of the calculation of such 
matrix elements has been described before for the case of 
Cs 3, 12, 1^]. However, we also find that term d, which 
is missing entirely from third-order calculation is rela- 
tively large for 7s — 5p and 7s — 7p matrix elements. It 
is rather puzzling since we observe very good agreement 
of third-order and all-order results for these transitions. 
The possible explanation is accurate cancellation of the 
high-order terms for these transitions. We note that term 
d is not significant for the 6s — 5p transition, for which 
some conclusion of accuracy may be drawn from the com- 
parison with similar Cs transition matrix elements (which 
are measured to quite good accuracy). Accurate mea- 
surements of these matrix elements would be very useful 
in developing further understanding of this issue. We 
note that even though this term is significant, it is not 
very large, on the order of 5%. 

In the case of nd — n'p transitions, the term c (or BO 
term for third-order calculation) dominates. The rela- 
tive contribution of the correlation correction varies with 
the transition. For 5c? — bp and 6c? — 6p transitions, the 
DHF approximation gives a very poor result; the sign of 
the matrix element changes when correlation is added for 
6c? — 6p matrix elements and the DHF values for 5c? — bp 
matrix elements are about 7 times too low. As expected, 
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TABLE III: Contributions to third-order and all-order electric-dipole reduced matrix elements (a.u.) for selected transitions 
in Rb. 





5s — 5pi 10 


5s — Tpi/2 


7s — 5pi/2 


7s — 7pi 10 5d-}/o 


— 5pi/2 


5d-j /9 — 7pi 10 


fSd-iio — 6pi 10 


6d.i /o — 4 fc 10 










Third order 










DHF 


4.819 


0.149 


0.981 


17.612 


0.244 


13.639 


-0.254 


6.109 


RPA 


-0.213 


-0.060 


0.015 


-0.024 


0.101 


0.030 


0.045 


0.002 


BO 


-0.419 


0.040 


-0.038 


-0.899 


0.888 


-3.048 


1.589 


3.041 


SR 


0.027 


0.009 


-0.004 


0.002 


-0.013 


-0.004 


-0.005 


0.001 


Norm 


-0.033 


-0.001 


-0.002 


-0.012 


-0.001 


-0.016 


0.000 


-0.004 


Total 


4.181 


0.130 


0.952 


16.679 


1.220 


10.600 


1.375 


9.150 










All order 










DHF 


4.819 


0.142 


0.981 


17.612 


0.244 


13.639 


-0.254 


6.109 


Term a 


-0.234 


-0.064 


0.014 


-0.025 


0.102 


0.036 


0.047 


0.005 


Term c 


-0.375 


0.009 


-0.077 


-0.648 


1.575 


-3.372 


2.752 


4.393 


Term d 


0.039 


0.017 


0.056 


0.149 


-0.197 


-0.024 


-0.439 


-0.089 


Other 


0.039 


0.012 


-0.006 


0.005 


-0.025 


-0.004 


-0.009 


0.002 


Norm 


-0.068 


-0.002 


-0.014 


-0.248 


-0.083 


-0.506 


-0.109 


-0.483 


Total 


4.220 


0.115 


0.954 


16.844 


1.616 


9.768 


1.989 


9.938 



we find very large, 25% — 30%, differences between tliird- 
order and all-order calculations for tfiese transitions as 
the all-order calculation includes correlation more com- 
pletely. The term c can be corrected by including triple 
excitations or use of the semi-empirical scaling described 
in Refs. 0, 0, ^3 ■ To check the validity of the scaling 
approach, we applied such a method to the 5c? — 5p ma- 
trix elements and compared the resulting ratio with the 
experimental measurement |14|. 



TABLE IV: The ratio i? of the 5ci3/2 — 5p3/2 and 5^3/2 — 5pi/2 
electric-dipole matrix elements divided by the corresponding 
values of the Ci reduced matrix elements, where Ci is the 
normalized spherical harmonic. Comparison of the present 
results calculated in different approximations with theoretical 
and experimental values from Ref. jlj. The corresponding 
5ci — 5p matrix elements, divided by the corresponding values 
of the Ci reduced matrix elements are also given. 



J- J' 


3/2 - 5/2 


5pj - 5dj' 
3/2-3/2 


1/2-3/2 


R 


DHF 


0.3182 


0.3034 


0.2115 


1.434 


Third order 


1.1756 


1.1758 


1.0562 


1.113 


All order 


1.5063 


1.5248 


1.3995 


1.089 


Scaled 


1.2801 


1.2882 


1.1707 


1.100 


MBPT [141 


1.0238 


1.0216 


0.9000 


1.135 


Expt. [14] 








1.068(8) 



In Table Hvl we give our values of the ratio R 



R = 



(5p3/2||-P||5d3/2) (5pi/2||Ci 115^3/2) 
(5pi/2||L'||5d3/2) (5p3/2||Ci 115^3/2) 



of the 5p3/2 — 5^3/2 and 5pi/2 — electric-dipole ma- 
trix elements divided by the corresponding values of the 
Ci reduced matrix elements, where Ci is the normal- 
ized spherical harmonic p^ . The present results are 
compared with theoretical and experimental values from 



Ref. |l4l| . The corresponding 5d — 5p matrix elements, 
divided by the corresponding values of the Ci reduced 
matrix elements are also given. The theoretical value 
of Ref. H was obtained using relativistic third-order 
many-body theory some higher-order Brueckner- 
orbital terms were also included, which accounts for the 
difference with our third-order result. As noted above, 
the largest contribution to the 5p — 5d matrix elements 
calculated using the all-order method comes from the 
term c containing single-valence excitation coefficients 
Pmv We scale these excitation coefficients pmv with the 
ratio of the experimental and correspondin g th eoretical 
correlation energies as described in Refs. '^,t5|. We list 
the results obtained using such scaling in a row labeled 
"Scaled" . We find very significant differences between all 
high-precision calculations, which indicates that the ac- 
curacy of the all-order calculation is around 10-20% for 
these matrix elements. However, our all-order result for 
the ratio R is in much better agreement with experiment 
01 than the third-order value and theoretical calculation 
of Ref. Q. 



III. LIFETIMES 

We use the resulting all-order matrix elements to cal- 
culate the lifetimes of the 6s, 7s, 8s, 6p, 7p, 8p, 4d, 5d, 
and 6d levels in Rb for the comparison with experiment 
and providing benchmark values for these lifetimes. 

The Einstein A-coefficients A.^w ^3 Eire calculated us- 
ing the formula 



2.02613 X lO^s KwllDllw)!^ 



A3 



2i. 



1 



(3) 



where (wHZJUw) is the reduced electric-dipole matrix el- 
ement for the transition between states v and w and A 
is corresponding wavelength in nm. The lifetime of the 
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state V is calculated as 



and the denominator of Eq. Q gives the radiative width 
of the level v. 



TABLE V: Transition energies (E/hc) in cm^^, matrix el- 
ements (in a.u.) and corresponding contributions to 7p3/2 
radiative width (in MHz). 



Transition (vw) 


{v\\D\\w) 




4 




7P3/2 - 7s 


23.587 


1559 


1.068 


1.068 


7P3/2 - 6s 


1.540 


7736 


0.556 


1.624 


7P3/2 - 5s 


0.202 


27870 


0.447 


2.071 


7^3/2 - 5^3/2 


4.242 


2169 


0.093 


2.164 


7p3/2 - 4(^3/2 


0.470 


8515 


0.069 


2.234 


7p3/2 - 5^5/2 


12.798 


2169 


0.847 


3.081 


7P3/2 - 4^5/2 


1.411 


8515 


0.623 


3.703 



Several electric-dipole (El) transitions contribute to 
the lifetime of each of the levels considered here. The 
simplest case is that of the 6s level, where only 6s — 5pi /2 
and 6s — 5p3/2 transitions need to be included. To calcu- 
late the lifetime of the 8^3/2 state, we need to include ten 
El transitions and in the case of the 6^3/2 level seven El 
transitions are allowed, including the 6^272" 4/5/2 transi- 
tion. The experimental energies from [ly are used in the 
lifetime calculation. We illustrate the importance of the 
different channels contributing to the lifetime of the 7p3/2 
state in Table EI where we give the coefficients Ay^ for 
each transition together with the transition energies from 
[Tll | and corresponding matrix elements from Tables HITll 
The accumulated sum A.^w is listed in the last col- 
umn. Wc find the contribution from the 7p3/2 — 5^^5/2 
transition to be nearly as large as the contribution from 
the 7p3/2 — 7s transition. All seven contributions need to 
be included in an accurate calculation. The smallest con- 
tribution to Ayw comes from 7p3/2 ~ 4^3/2 transition 
and is around 2%. 

We list the contributions to 6^3/2 and 60/5/2 lifetimes 
as well as the corresponding matrix elements and tran- 
sition energies in Table IVTI The dominant contributions 
to 6^3/2 and 6^5/2 lifetimes come from the 6(^3/2 — 5pi/2 
and 6c?5/2 — 5p3/2 transitions, respectively. The next 
largest contribution is 14% for the 6^3/2 lifetime and 7% 
for the 6(i5/2 lifetime. We find the contributions from 
6*^3/2 ~ 4/5/2 and 6^5/2 — 4/7/2 transitions to be around 
7% for the corresponding level. 

The results for the ns, np and nd lifetimes obtained 
using all-order matrix elements are compared with ex- 
perimental values from Refs. 0,0, in Table Ivnl We 
omit the effect of the blackbody radiation in our calcu- 
lations. In Ref. 101, it was estimated to be small for the 
levels considered here with the exception of the 8p levels. 
It is still far below experimental uncertainty of the 8p 
level lifetime. The results agree with experiment within 



TABLE VI: Transition energies E/hc in cm^^, matrix ele- 
ments (in a.u.) and the corresponding contributions to 60(3/2 
and 6(^5/2 radiative widths (in MHz). 



Transition [vw) { 


v\\D\\w) 


3 Euyj 


Ayw y ^VW 


6^3/2 - 5pi/2 


1.180 


16108 


2.948 2.948 


6^3/2 - 6pi/2 


1.989 


4972 


0.246 3.194 


6^3/2 - 7pi/2 


31.422 


852 


0.309 3.504 


0"3/2 JP3/2 




LOO 1 \J 




6^3/2 - 6p3/2 


1.012 


4894 


0.061 4.195 


6^3/2 - 7p3/2 


14.161 


817 


0.055 4.250 


6^3/2 - 4/5/2 


9.938 


1895 


0.341 4.591 


6d5/2 - 5p3/2 


1.658 


15872 


3.712 3.712 


6^5/2 - 6p3/2 


2.974 


4896 


0.351 4.063 


6^5/2 - 7p3/2 


42.481 


819 


0.335 4.398 


6^5/2 ~ 4/5/2 


2.642 


1897 


0.016 4.414 


6^5/2 - 4/7/2 


11.813 


1897 


0.322 4.736 


TABLE VII: Radiative widths (in MHz) 


and lifetimes (in ns) 


for ns, np, and nd states in Rb. 


Experimental values are taken 


from the compilation 


in Ref. 15] unless noted otherwise. The 


lifetimes obtained using ab initio all-order matrix elements 


and scaled all-order values are 


listed in 


columns Tv and t^"^, 


respectively. 








State ^ ] 




' V 


cxpt 

' V 


6s 21.761 


46.0 


45.4 


46(5) 


7s 11.318 


88.4 


88.3 


88(6) 








91(11) 


8s 6.201 


161.3 


161.8 


161(3) [4] 








154(7) 








153(8) 


6pi/2 7.729 


129 


123 


125(4) 








131(5) 


7pi/2 3.394 


295 


280 


272(15) 


8pi/2 1.899 


527 


508 




6p3/2 8.461 


118 


113 


112(3) 


7p3/2 3.703 


270 


258 


246(10) 








233(10) 


8p3/2 2.056 


486 


471 


400(80) 


4d3/2 11.478 


87 


83.5 


86(6) 


5d3/2 5.193 


193 


243 


205(40) 


6d3/2 4.591 


218 


263 


256(4) [4] 






298(8) [5j 


4d5/2 10.674 


94 


90 


94(6) 


5d5/2 5.340 


187 


235 


230(23) 


6d5/2 4.736 


211 


252 


249(5) [4] 


the experimental precision for 6s, 7s, 


8s, 6p, and Ad lev- 


els. The values of the 7p, i 


]p, 5d, and 6d lifetimes are 



in good agreement with experiment. The lower accuracy 
of the calculations is expected for these levels owing to 
the lower accuracy of the np — n'd matrix elements with 
comparison to most of the np — n' s ones due to large cor- 
relation corrections for these matrix elements. It should 
be noted that the experimental uncertainties of 7p, Bp, 
and Sd lifetimes are rather large (up to 20%) and two 
recent experiments for the 61^3/2 level give substantially 
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different values: 256(4) ns and 298(8) ns. Our result 
for 6(i3/2 lifetime (263 ns) is very close to the first mea- 
surement, however, large uncertainties in the correlation 
contributions to the relevant transitions do not allow to 
rule out 10% larger value of this lifetime. 

The accuracy of the theoretical calculations for these 
levels may be improved by including some triple- 
excitation contributions into the corresponding wave 
functions. 

We estimate some omitted higher-order contributions 
using the scaling of the single excitation coefficients Pmv 
described above. The contributions containing these 
single-excitation coefficients are dominant for all of the 
transitions needed for the calculation of lifetimes in Ta- 
ble IVIII The lifetimes obtained using scaled all-order 
matrix elements are listed in column r^'^. We find that 
the scaled data agree with experimental values within the 
experimental uncertainty for all levels with exception of 
7p3/2 and 6^3/2 , where the theoretical values are just 
outside the experimental uncertainties. 



-0.5 









5s-5p^^2 







0.056 0.057 0.058 0.059 0.06 

CO (a.U.) 



FIG. 1: Dynamic polarizability a{uj) for the ground state of 
Rb in a.U. 



IV. POLARIZABILITIES 

The valence part of the ac polarizability of an alkali- 
metal atom in ns state can be calculated using the for- 
mula 

. ^ 1 / (En'py, - Ens){n'pi/2\\D\\nsy 

(i?„-^3,,-i?„,)2-c.2 j' ^> 

where D is the dipole operator. In this formula, oj is 
assumed to be at least several linewidths off resonance 
with the corresponding transition. 

The core contribution to the polarizability, calculated 
in the DHF approximation is found to be small for Rb 
(9.3 Oq) and is weakly dependent on uj in the frequency 
range considered here. The static value for the polariz- 
ability of Rb+ calculated in the random-phase approx- 
imation (RPA) ^6, 17] is 9.1 an, close to the value of 
9.0 ag obtained by Johansson jl^ from analysis of the 
observed term values of nonpenetrating Rydberg states. 
The accuracy of the RPA approximation for the core po- 
larizability is estimated to be 5 % in Ref. Q . We use the 
RPA value for the core polarizability of Rb"*" as a base- 
line, and adjust it to account for valence electron (using 
Eq.© with n' — 2,3,4) and the frequency dependence 
by using DHF calculations. The RPA and DHF values 
differ by only 2 %. The correction to the core polariz- 
ability owing to the presence of the valence electron is 
very small, it is only — 0.3aQ for 5s state in DHF approx- 
imation. 

First, we describe the calculation of the ground state 
polarizability, a^si^)- The expression © converges 
rapidly with n' so the contribution with n' = 5 is domi- 
nant. We use experimental 5s — 5p matrix elements from 



0, all-order matrix elements from Ref. and experi- 
mental energies from to evaluate the expression of 
Eq. © with n' = 5, 6, 7, 8. The contribution to the 
ground state polarizability from states with n' > 8 is very 
small, 0.2aQ in DHF approximation. We plot the ground 
state frequency-dependent polarizability in the vicinity 
of the 5s — 5pi/2 and 5s — 5p3/2 resonances in Fig.^ The 
behavior within the few linewidths from the resonances 
is not shown and the exact placement of the resonances 
is indicated by the vertical dashed lines. As we see from 
Fig. n a5s(w) changes sign between two resonances. We 
determined that a^s{^) = at A^ac = 790.032(8) nm 
^2|, where the uncertainty results from the uncertainty 
of the polarizability calculation. This crossing point is 
of interest for the optical lattice experiments with Rb as 
the atoms will no longer be trapped at this wavelength. 

We give the breakdown of the different contributions 
to the ground state polarizability for one particular fre- 
quency, uo = 0.04298 a.U., corresponding to A = 1.06/im 
, in Table IVIIII The comparison with other theory 
|20j | and experiment |2lj is also given. Our value is in 
good agreement with the result from ^2^. The discrep- 
ancy is due to our use of more accurate values for the 
electric-dipole matrix elements. Our result is just out- 
side the range of uncertainty of the value a = 769 ±61 Oq 
inferred by Bonin and Kadar-Kallen from an atomic de- 
ffection experiment [2l| . 

Next, we describe the calculation of the polarizabilities 
of the Rydberg states. We calculate the polarizabilities 
of the ns states of Rb with n = 8 ... 19 using the DHF 
approximation; i.e. using DHF values for both energies 
and matrix elements in Eq. ((S)). The summation over n' 
in Eq. is truncated at n' = 23. The DHF calculations 
are done on the non-linear grid of the form 

r{i) = ro(e(^-i)'' - 1). (6) 

To calculate matrix elements and energies of the first few 
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TABLE VIII: Contributions to ground state dynamic polar- 
izability of Rb in for uj = 0.04298 a.u. 



Contribution 


Value 


1 r 

Tl — 




n = 6 


i.4o 


n — / 


U.iO 


n' = 8 


0.04 


a„(n' = 5..8) 


684.52(72) 


Q„(n' > 8) 


0.17(8) 


a„c(n' = 2. .4) 


-0.26(13) 


etc 


9.10(45) 


Final a 


693.5(9) 


Ref. [22] 


711.4 


Expt. [21J 


769(61) 


TABLE IX: Dynamic polarizabilities a„s{i^) (in units of ajj) 


for Rb, cj=0.0576645 a.u., i.e. A = 790 nm. 




n Qdhf Q 


n Qdhf 


8 -304 -295" 


14 -286 


9 -292 


15 -285 


10 -289 


16 -284 


11 -288 


17 -282 


12 -287 


18 -280 


13 -287 


19 -277 



"High-accuracy value obtained using experimental energies and 
all-order matrix elements for the dominant terms with n' = 7, 8. 



ered here. The polarizabilities of the ground and Ry- 
dberg states can be matched at the point between two 
resonances where the ground state polarizability is small 
and negative. The exact matching point for the 15s state 
is Xvac = 790.14(2) nm 2]. The matching of the po- 




0.056 0.057 0.058 0.059 
co(a.u.) 

FIG. 2: Dynamic polarizability a{ijj) for the ns states of Rb 

in Oq. 



larizabilities allows to minimize motional heating in the 
quantum computation scheme with neutral atoms. 



excited states 500 grid points are sufficient. This type of 
grid provides a very dense grid near the origin. Thus 
we needed to confirm that there were sufficient points at 
large R to support the Rydberg states. The parameters 
ro and h were chosen to ensure sufhcient number of points 
to calculate {ns\\D\\n'p) matrix elements for high values 
of n and n' . However, we found that the accurate DHF 
values of the matrix elements even between high Ryd- 
berg states can be obtained with relatively small number 
of grid points. For example, the values of the 20s — 20p3/2 
matrix element obtained with 500 and 40000 point grids 
(corresponding to approximated 10 and 200 points in the 
relevant grid section) differ by only 0.05%. The results of 
the DHF calculation for tj=0. 0576645 a.u. corresponding 
to A = 790 nm are given in Table IIXI We find that the 
polarizability values do not change substantially with n 
for n > 8. We conduct a separate high-precision calcu- 
lation of the 8s polarizability using the all-order matrix 
elements from Table We find DHF result to be in good 
agreement with the high-precision value. 

We plot the frequency-dependent polarizabilities of 8s, 
15s, and 20s states in Fig. |21 As expected, the polariz- 
abilities depend weakly on uj for the frequencies consid- 



V. CONCLUSION 

We have conducted a systematic study of the ns — n'p 
and nd — n'p electric-dipole matrix elements in rubid- 
ium using relativistic all-order method. An investigation 
of the accuracy of these matrix elements was performed. 
The resulting matrix elements were used to calculate life- 
times of the ns and np levels with n = 6, 7, 8 and nd levels 
with n = 4, 5, 6. The lifetime values were found to be in 
good agreement with experiment. The dynamic polar- 
izabilities of the ns Rb states, which are of interest for 
the optimization of quantum computation scheme with 
neutral atoms mediated by the conditional excitations to 
Rydberg states, were also calculated. 
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